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Abstract 

Let H(h) = —h 2 d 2 /dx 2 + V(x) be a Schrodinger operator on the real line, 
W(x) be a bounded observable depending only on the coordinate and A; be a 
fixed integer. Suppose that an energy level E intersects the potential V{x) in 
exactly two turning points and lies below Voo = liminf|, I .|^ 00 V(x). We consider 
the semiclassical limit n — > oo, h = h n — > and E n = E where E n is the 
nth eigen-energy of H(K). An asymptotic formula for (n|W(a;)|n + k), the non- 
diagonal matrix elements of W(x) in the eigenbasis of H(h), has been known 
in the theoretical physics for a long time. Here it is proved in a mathematically 
rigorous manner. 

Keywords: semiclassical limit, non-diagonal matrix elements, WKB method 

1 Introduction 

In the quantum mechanics the matrix elements of an observable occur in various 
situations. Let us mention few of them. They measure transition probabilities between 
two states and the coefficients in the stationary perturbation theory are expressed in 
terms of the matrix elements of the perturbation. The distribution of matrix elements 
is of interest for quantum systems stemming from classically chaotic systems, see 
for example jH [6] and references in the latter paper. Our immediate motivation to 
study the matrix elements was the quantum version of the Kolmogorov-Arnold-Moser 
method pQ, [8]. One of the assumptions under which this method is applicable is that 
a time-dependent perturbation of a quantum system must be sufficiently small with 
respect to certain norm which is also expressed in terms of matrix elements. 

One may hope to obtain at least a qualitative information about the behavior of 
matrix elements when considering the semiclassical limit. In fact this idea goes back to 
the very origins of the quantum mechanics. A semiclassical formula for non-diagonal 
matrix elements in the one-dimensional case has been suggested already a long time 
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ago [12]. In [H] one can find another derivation, also on the level of rigor usual in the 
theoretical physics, for absolute values of the non-diagonal matrix elements. 

Despite of the ancient history rigorous mathematical results have been published 
essentially more recently. Moreover, they cover only some particular cases even though 
the technical tools necessary for the derivation may be at hand nowadays. One usually 
assumes that the corresponding classical system is either ergodic [5] , [B] or completely 
integrable [TH], [2], [15], [7J. The semiclassical limit of diagonal matrix elements is now 
treated in detail [5]. In the case of multi-dimensional completely integrable systems a 
formula for non-diagonal matrix elements was proved in [19], [15], [7J, see also [16] for 
some generalizations. The one-dimensional case seems to be rather particular. In [H] 
one can find a derivation of the semiclassical formula for pseudo-differential operators 
in one variable such that the Weyl symbol of the Hamiltonian is a real polynomial on 
the phase space while imposing an additional assumption on the discreteness of the 
operator spectrum. 

The present paper aims to provide a mathematically rigorous verification of the 
semiclassical limit of non-diagonal matrix elements for Schrodinger operators on the 
real line. We prove the formula under mild assumptions on the potential. In addition, 
we take care about identifying the quantum number coming from the Bohr-Sommerfeld 
quantization condition with the index determined by the natural enumeration of eigen- 
values in ascending order. Our approach relies on a transparent application of some 
well established tools in the spectral and semiclassical analysis. So we briefly recall the 
corresponding results while adjusting their formulation to our purposes. On the other 
hand, the chosen method restrict us to considering observables which depend on the 
coordinate only. This particular case was sufficient for the applications we originally 
had in mind, as mentioned above. 



Let us now formulate precisely in what sense the semiclassical limit is understood. 

Set 



We consider a fixed energy E and an observable W = W(x) depending only on the 
coordinate x. The assumptions are as follows. 

We suppose that V(x) is bounded from below and three times continuously differ- 
entiable, W(x) is bounded and continuously differentiable, 



We assume that at the energy E there are exactly two regular turning points, i.e., 
V~\E) = {£_,£+}, x_ < x+, and V'(x±) ^ 0. Set 




(1) 




(2) 




(3) 



In addition we introduce an assumption making it possible to apply the WKB approx- 
imation, namely we assume that 



where a is a positive number chosen so that f(x) > 5 > for \x\ > a. Notice that 

1 d 2 / 1 \ b(V') 2 — 4(V — E)V" 



fV*dx 2 V/ 1/4 / ~ 16(V-E) 5 / 2 

It may be convenient to replace condition (J4j) by two simpler conditions, 

\V'\ 2 f \V"\ 

dx < oo, / — — -di<oo. (5) 



-a,a] ~ E) 5 / 2 ' J m . a , a] (V - Ef/1 

The part of the spectrum of H(h) lying below is known to be formed exclusively 
of simple isolated eigenvalues. We fix the phase of an eigenfunction ip n corresponding 
to an eigenvalue E n < by requiring ip n to be positive on a neighborhood of +oo. 
Moreover, there exists a strictly decreasing sequence of positive numbers tending to 0, 
{fcn}%Ln i an d a constant ho > such that for h e ]0, Hq], E belongs to the spectrum 
of H(h) if and only if h = H n and in that case E = E n is the nth eigenvalue of H{K) 
provided the enumeration of eigenvalues starts from the index n = 0. 

Under these assumptions we claim that if k e Z is fixed, n — > oo, K = h~ n — > 0, 
with E = E n , then 

(n\W{x)\n + k) -> i / W(q(t))e ikuJt dt (6) 

Jo 

where (q(t),p(t)), t G [0,T], zs i/ie classical trajectory in the phase space at the energy 
E and with the initial point chosen so that the kinetic energy vanishes, i.e., p(0) = 0, 
and g(0) coincides the right turning point x + . Furthermore, T > is the period of the 
classical motion and u = 2-k/T is the frequency. 

Remark. If the phase of the wave function ip n was chosen so that ip n was positive on a 
neighborhood of — oo then formula ([6]) would be again true with (g(0),p(0)) = 0). 

As already said, we have confined ourselves to observables depending only on the 
coordinate because our method of proof is based on the WKB approximation. One 
naturally expects, however, that for any smooth bounded classical observable A(q,p), 

(n\A\n + k) - ^ T A(g(t),p(t))e^dt 

where A is a suitable quantization of A. We have already mentioned that this result 
is actually proved in [TJ] in the case when the potential V(x) is a polynomial. 

Let us rewrite the RHS in formula (E]). The equation of the classical trajectory in 
the phase space reads p 2 + V(x) = E and its period equals 

x+ dr 

(7) 



\[E^V{ 



x 



For x E [x- , x+] set 



i i 1 I'* d V /m 
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Then t(x + ) = 0, t(x_) = T/2, g(r(x)) = x, and 



/ W(q(t))e ikujt dt = [ 



■x+ 



W{x) 



2itk 



t(x) J dx. 



VE - V(x) 



cos 



T 




The paper is organized as follows. In Sections 2 through 4 we recall some prelim- 
inaries that we need for the proof of the formula. Section 2 is devoted to the basic 
spectral properties of the Schrodinger operator, Section 3 is concerned with the Weyl 
asymptotic formula and some basic facts about the WKB approximation are sum- 
marized in Section 4. By counting the zeroes of wave functions we show in Section 5 
that the quantum number coming from the Bohr-Sommerfeld quantization condition 
equals the index of the corresponding eigenvalue. The semiclassical formula is then 
proved in Section 6. 

2 Properties of the spectrum lying below 

Here we briefly recall two well known properties of Schrodinger operators. In the mono- 
graphs they are usually formulated and derived for potentials diverging at infinity. We 
just wish to point up that the same assertions apply also for more general potentials 
provided one takes care only about the part of the spectrum lying below V^. The 
corresponding proofs can be taken almost literally from the cited monographs. 

In this section (and only in it) the Planck constant is not relevant and so we set it 
equal to 1 and consider the Hamiltonian 



The following theorem is in fact widely used. We recall it in a form which is a direct 
modification of Theorem XIII. 16 in [T7]. Its proof is based on the min-max principle 
and is applicable in any dimension of the underlying Euclidean space. Moreover, the 
differentiability of V(x) is not required. 

Theorem 1. Let V be a measurable function in R™ which is bounded from below. De- 
fine H = — A + V as the sum of quadratic forms in L 2 (R™, d n x). Then the lower edge of 
the essential spectrum of H , if any, is greater than or equal to = lim iniM-^ V(x). 

Let us note that in the one-dimensional case and provided the potential is contin- 
uous Theorem [T] also follows from a well known estimate on the number of negative 
eigenvalues. 

Here and everywhere in what follows, if A is a self-adjoint operator then P(A; •) 
designates the associated projector- valued measure, and for K e M we denote 



H 



dx 2 



+ V(x) inL 2 (M,dx). 



N(A,K) =TimkP(A;}-oo,K[). 



Further, for a real- valued function W(x) we set 



W-{x) = max{0, -W(x)}. 
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It holds (see, for example, Theorem 5.3 in [3]) 

N(H,0) < 1+ / \x\V^{x)dx. 
Jr 

In particular, if V(x) is continuous and bounded from below then for any c < the 
function (V — c) _(x) has a compact support and, by this estimate, N(H, c) < oo. This 
again implies that the lower edge of the essential spectrum of H is greater than or 
equal to V^. 

The next property is specific for the one-dimensional case. The potential V(x) is 
supposed to be continuous and bounded from below. 

As is well known from the theory of ordinary differential equations, for E < V^, any 
nontrivial solution of the Schrodinger equation either grows at least exponentially or 
decays at least exponentially at +oo (see, for example, Corollary 1 in [HI Section II]). 
The latter solution is called recessive at +oo and is unique up to a multiplicative 
constant. Of course, an analogous assertion is also true for — oo. It immediately follows 
that all eigenvalues of the Hamiltonian H lying below are simple. Moreover, in 
virtue of Theorem [H they have no accumulation points below V^. Consequently, the 
eigenvalues of H below can be arranged into a strictly increasing sequence, empty 
or finite or infinite, 

E Q < E x < E 2 < . . . <Voo. 

The following theorem is a straightforward modification of Theorem 3.5 in [31 
Chapter II]. 

Theorem 2. The number of zeroes of the mth eigenfunction of H corresponding to 
the eigenvalue E m < is exactly equal to m. 

3 The Weyl asymptotic formula 

In this section we aim to recall the Weyl asymptotic formula generalized to Schrodinger 
operators. It can be derived from the Gutzwiller trace formula [TU] which was rigor- 
ously proved in |4| under the assumption that the potential is positive and infinitely 
differentiable. In [18] there is given a short review of the history and the Weyl asymp- 
totic formula is recalled even under stricter assumptions which among others mean 
that the potential does not grow faster than polynomially. A weaker version of the 
formula is also stated in [T71 Theorem XIII. 79] but only for compactly supported 
potentials. 

Here we wish to point out that the proof of Theorem XIII. 79 in [17] can be extended 
in a straightforward manner and thus the Weyl asymptotic formula can be derived just 
under the assumption that the potential is semi-bounded and continuous. We restrict 
ourselves, however, to the one-dimensional case only. In addition, this approach is 
quite simple as it is based merely on an application of the min-max principle and the 
Dirichlet- Neumann bracketing. On the other hand, if compared to the result based on 
the trace formula, as presented in [18], the control of the error term is essentially worse; 
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it is known to be of order 0(1) while the present method only yields the asymptotic 
behavior of the type o{fT l ). 

From now on, the Planck constant is again relevant. This means that the discussion 
concerns the Hamiltonian H(h) introduced in ([1]). Since what follows is nothing but 
a slight modification of known results we just indicate the basic steps. 

First let us recall a definition from [T7[ XIII. 15] making it possible to compare 
self-adjoint operators defined in different Hilbert spaces. The symbol Q(A) stands for 
the form domain of A. If if? G Q{A) then the scalar product (if?, Aif>) is automatically 
understood in the form sense. 

Definition. Let Hi C H be a closed subspace, let A be a semi-bounded self-adjoint 
operator in H and let B be a semi-bounded self-adjoint operator in H\. We shall write 
A < B if and only if it holds 

(i) Q(A) d Q(B), 

(ii) Wif?eQ(B), (if;,Aif;) < (if;,Bif;). 

With the aid of the min-max principle one can show [TTl XIII. 15] that if A < B 
then 

(i) VK G R, rank P(A; } - oo, K[) > rankP(P; ] - oo, K[), 

(ii) VK G R, rankP(A; ] - oo, K]) > rankP(P; ] - oo, K]). 

The following lemma is analogous to Proposition 2 in [TTl XIII. 15] in the one- 
dimensional case and its proof is based on rather elementary explicit computations of 
the eigenvalues for the involved operators. 

Lemma 1. Let I = [a, b] be a compact interval. Let us introduce Ho, Hn and Hm as 
self-adjoint operators in L 2 (I,dx) such that all of them act as the differential operator 
—h 2 d 2 /dx 2 and whose domain is respectively determined by the Dirichlet, Neumann 
and mixed boundary conditions. Then for all K > it holds 

-1 < rank P(H;] - oo,K[) -VK < rank P(H; ]- oo, K\) -VK < 1, 

Tin Tin 

where H is any of the operators H D , H N , H M , and t = b — a is the length of the 
interval. 

The following lemma coincides with Proposition 4 in [TTl XIII. 15] in the one- 
dimensional case. 

Lemma 2. Let — oo < a < b < c < +oo and let H be a self-adjoint operator in 
L 2 ([a,c],dx) which acts as the differential operator — d 2 /dx 2 with either the Dirichlet 
or the Neumann boundary condition imposed at each of the points a and c (mixed 
boundary conditions are admitted). Let H$ and H$ be the self-adjoint operators in 
L 2 ([a, b], dx) also acting as — d 2 /dx 2 and with the domain being determined by the same 
boundary condition at the point a as imposed in the case of the operator H and by the 
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Dirichlet or Neumann boundary condition at the point b, respectively. Analogously one 
introduces the self-adjoint operators Hp and in L 2 ([b, c], dx). Then it holds 

H$ © <H< 4 1 ' © H%\ 

First let us state the Weyl asymptotic formula for a finite interval. It can be prover 
in a way very close to the proof of Theorem XIII. 79 in [T7]. So we do not reproduce 
the proof but let us note that it is based on a limit procedure when the interval is 
split into N sub intervals of equal length with N tending to oo. In the course of the 
proof one uses Lemma Q] and [2J the additivity of the numbers N(A,K), i.e., 

N{A 1 © A 2 © . . . © A N , K) = N(At, K) + N(A 2 , K) + . . . + N(A N , K), 

and the fact that the integral on the RHS of (J9]) exists in the Riemann sense. 

Theorem 3. Let — oo < a < b < +00, V G C([a, b]), and let 

be a self-adjoint operator in L 2 ([a, b], dx) with either the Dirichlet or Neumann bound- 
ary condition imposed at each of the boundary points a and b (mixed boundary condi- 
tions are admitted). Then for all K G R, 

1 f b 

lim HN(Hf(h),K) = - / ^f(V -K)_(x) dx. (9) 

H->0+ ' 7T J a 

Finally let us proceed to the case of the Hamiltonian H(h). 

Theorem 4. Let V G C(R) be a real-valued function which is bounded from below. 
Then for all K < Voo it holds true that 



lim hN{H{h),K) = — Vol z -oo,K[)) = - [ y/(V - K)-(x) dx 



(10) 



where Jif{x,p) = p 2 + V(x) and Vo\z(X) designates the Lebesgue measure of a mea- 
surable set X in the phase space. 

Proof. If K < Voo then the support of (V — K)_ is compact. Suppose that 
supp(V — K)_ C [a, b], —00 < a < b < +00. Set 

d 2 

HAh) = -h 2 — - (V-K)Jx) inL 2 (M,dx) 
dx 1 

and 

d 2 

H 2 (h) = -h 2 — + V(x) - K in L 2 ([a,b},dx) 

with the Dirichlet boundary condition imposed at the points a and b. Observe that 
— (V - K)-(x) < V(x) - K on R and so Q(H(h) - K) C Q(H x {h)). Furthermore, 
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L 2 ([a,b],dx) can be naturally regarded as a subspace in L 2 (IR, dx). If ip G Q(H 2 (h)) 
then ip defined by ip(x) = ip(x) for x G [a, b], ip(x) = for x G R \ [a, 6], belongs to 
Q(H(h) — if) ("0 is an absolutely continuous function). This implies that Q(H2(h)) C 
Q(#(fr) - if). We have find that H^h) < H(H) -K< H 2 (h). Hence 

N(H 2 (h),0) < N{H{h),K) < iV(tfi(fc),0). 

Formula ( TTU]) for compactly supported potentials is stated in jTTl Theorem XIII. 79]. 
Hence it holds 

lim HNiHtih)^) = - I ^{V -K)^{x) dx, 
and from Theorem [3] we know that 

lim HN(H 2 (h),0) = - [ y/(V-K)_(x) dx = - [ y/(V-K)_(x) dx. 

h ^°+ ' K J a 7T J R 

Formula ( jTUl) for a general potential then follows by bracketing. □ 

For our purposes the following immediate corollary of Theorem H] will be sufficient. 
Suppose that V(x) is continuously differentiable and an interval ]a, b[, a < b < Voo, 
contains at least one regular value of the classical Hamiltonian Jf?(x,p), i.e., there 
exists A e]a,b[ satisfying J^ _1 ({A}) ^ and V(x) = X implies V'(x) ^ 0. Then the 
number of eigenvalues of H(h) in the interval ]a, b[ tends to infinity as h — > 0+. 



4 The WKB method for one and two turning points 

Here we summarize some basic facts about the WKB approximation, also called 
Liouville- Green approximation, that we need for the proof of the formula in Section [6j 
At the same time we introduce the necessary notation. We stick to the presentation 
given in the monograph [JS] whose distinguished feature is that it provides explicit 
bounds on the error terms. 

Let us first consider the situation with one turning point. Let ]a, &[C R be an 
interval, finite or infinite, xo G }a, b[, and f(x) be a real-valued function defined on 
}a,b[ such that f(x)/(x — x ) is positive and twice continuously differentiable (hence 
f(x ) = 0, f'(x ) > 0). For x G }a, b[ set 

2 
3 



C 3/2 = f* VWjto if^>^o, (Ha) 

\{-cf /2 = Hv^rndt ifx<x . (lib) 

Then ((x) is strictly monotone, ((x)/(x — Xo) is positive and twice continuously dif- 
ferentiable in ]a,b[, see Lemma 3.1 in [T^ Chapter 11]. 
Assume further that 



V / 7Wdt = oo (12) 
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and 



/ 

J]a. 



b[\Ui 



\f"\ f if') 2 

1 di<oo, / ^^dt<oo, (13) 



l/! 3/2 ' J]a M \U \f\ 5/2 



where U = [x — e,x + e] and e is any positive number such that a < x — e and 
x + e < b. 

Notice also that ^ 

C = (£) and C'(xo) = /'(x ) 1/3 . (14) 

Denote by £ the inverse function to £. Theorem 3.1 in [T3], Chapter 11, §3.3] can be 
rephrased as follows. 

Theorem 5. Under the above assumptions, the solution of the differential equation 

h 2 ^ = f(x)w (15) 

which is recessive as x tends to b exists on ]a, b[, is unique up to a multiplicative 
constant and equals 

1>(?) = (^ 1/4 (Ai(r 2 / 3 C) + e(h, x)) (16) 
with the error term satisfying the estimates 



\e(h,x)\<%(h- 2 ^0 h, 



de(h, x) 



dx 

where $>q(x), <&i(x) are certain continuous positive functions on R such that 



exp(-|x 3/2 ) 
const ^-7- as x ^ +00, 

2-1/4 

1 

const , 7- as x — > —00, 



\x 



1/4 



const exp^—-x 3 ^ 2 ^ as x — > +00, 
const as x — > —00. 



Let us now turn to the case when /(x) is given by (jHj) and so is defined on the 
entire real line. From now on the potential V satisfies all assumptions as formulated 
in the Introduction. In particular, it follows that the function 

V(x) — E 

is positive on R and belongs to C 2 (IR). (17) 



(x — X-)(x — x + ) 
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Moreover, there exists an open neighborhood of E, Ue =]E-, E + [, EL < E < E + , 
such that these assumptions apply for any A G Ue as well. 
For A G U E set 

7a = ^-\{X}) 

where Jif(x,p) = p 2 + V(x). Thus 7a is a closed curve in the phase space and the 
energy takes on it the value A. Let us further introduce the action integral, 



J(A) = f dxdp= [ pdx = 2 [ ( } y/X - V(x) dx 

J.je(x,p)<\ J ~ix Jx-{\) 



where x_(A) < x + (X) are the turning points at the energy A. Then 

T(X) = J\X) = / (19) 

is the period of the classical trajectory in the phase space. 

In the following theorem we summarize the result derived in [131 Chapter 13, §8.2]. 

Theorem 6. Under the assumptions on V formulated in the Introduction (in particu- 
lar, we assume that condition p7| ) is fulfilled as well as the convergence of the integrals 
in ^)) there exist a neighborhood Ue of E, h Q > 0, n G N and for every X e U E a 
sequence {ft n (A)}£L no , h > h no (X) > H no+ i(X) > h no+2 (X) > ... > 0, such that for 
h g]0, Hq[ the energy A is an eigenvalue of H(h) if and only if h = h n (X) for some 
n > n . Moreover, the sequence {^n(A)} asymptotically behaves like 

hniX)- 1 = (2n + I)ttJ(A)- 1 + 0(n- 1 ) (20) 

where the error term 0(n~ l ) decays in n uniformly with respect to X £ Ue. 

Remark. It is known that if V G C r (M), with r > 1, and E < Vac is a regular value of 
V(x) then the action integral J{X) defined in ( fl8l) is r times continuously differentiable 
on some neighborhood of E (see, for example, [18J). 

The verification of this assertion is quite elementary in the one-dimensional case 
and with two turning points at the energy E. For a sufficiently small neighborhood 
Ue = ]E-, E + [ the function V(x) is strictly decreasing on the interval [x- (E + ), X- 
and strictly increasing on [x+(E-),x+(E + )], with nowhere vanishing derivative. Let 
us write 



T(A) 




X+(A) \ dx 



a/A - V(x) 



T_(A)+T (A)+T + (A). 



Clearly, T (A) G C°°(U E ). Thus it is sufficient to verify that 71(A), T+(A) G C^^Ue) 
Let us focus only on the latter function. Set W+ = \ V\, ,„ . ,„ ) . Hence W+ 
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is r times continuously differentiable. After some elementary manipulations one can 
show that 



T+(A) 



f +(A) dx = 2 fx^W f - dt 

Jx + (E.) JX-V(x) v Jo V' 



x+ (e.) ^X - V(x) V Jo V'(W + (X(l-t^) + E^)) ' 



From the last expression it is obvious that T + (A) is r — 1 times continuously differen- 
tiable. 



5 Number of zeroes derived from the WKB method 

We need to show that if h = h m (X) and hence A is an eigenvalue of H (H), as claimed in 
Theorem El then A is exactly the mth eigenvalue of H(K). According to Theorem [21 the 
index of an eigenvalue lying below Voo equals the number of zeroes of the corresponding 
eigenfunction. Fortunately, the WKB approximation, as explained in [T3], is precise 
enough to control the number of zeroes. 

Let us recall some facts concerning the Airy functions. Let us denote by a n and b n 
the zeroes of the Airy functions Ai(x) and Bi(x), respectively, arranged in ascending 
order of the absolute value, i.e., . . . < b 3 < a 2 < b 2 < a\ < b x < 0. It is known that 

1\ _/ 1\\ 2/3 , /3 / 3\ _/ 3 nn2/3 



" (21) 
where ~b(x) = 0(x 1 ). 

First we again consider the situation with one turning point. Recall defining re- 
lations (II lap . (Illbl) for (. In the following theorem we summarize the results from 
§§ 6.1, 6.2 and 6.3 in pH Chapter 11]. 

Theorem 7. Under the same assumptions as in Theorem let w(x) be a nonzero 
solution of the differential equation 173]) on }a,b[ which is recessive as x tends to b 
(hence w(x) is unique up to a multiplicative constant). Then the set of zeroes ofw(x) 
in ]a, b[, denoted {z n } n >i and arranged in descending order, is at most countable. Any 
such a zero z fulfills £(z) < h 2 ^ 3 bi. Furthermore, for all sufficiently small h~ it is true 
that if C( a ) < fr 2 ^ 3 b n+ i then the nth zero, z n , does exist and obeys the estimate 

h 2/3 b n+1 < C(z n ) < h 2 ' z b n . 

Moreover, it holds 

\C{z n ) - 3 a n \ = 0(rT^)h 
where the symbol 0(n~ 1 ^) is uniform with respect to h. 

Remarks. From Theorem [7] it immediately follows that there are no zeroes in the 
interval [xo,6[. Furthermore, the number of zeroes of w(x) in any fixed nonempty 
subinterval ]c, d[ c]a, Xq[ tends to infinity as h — > 0+. 
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Now we come back to the case when f(x) is given by (JS]), with V{x) satisfying the 
assumptions from the Introduction. In particular, there are two turning points at the 
energy E, x_ and x + , and V(x) satisfies (IT7j) and (jSJ). Then for any a, x_ < a < x + , 
the function f(x) satisfies the assumptions of Theorem [7] with b = +00 and Xq being 
replaced by x + . Actually, condition (jSJ implies ffTB"]) and condition ffT2l is fulfilled 
automatically for < V^. Analogous arguments apply also for the other turning 
point x_. 

According to Theorem [6] there exist ho > and a sequence {h n }^L no , ho > h no > 
h no+ i > h no+2 > • • • > 0, such that for h G ]0, h [, E is an eigenvalue of H{h) if and 
only h = h n for some n > n . Let ip n (x) be an eigenfunction of H(h n ) corresponding 
to the eigenvalue E. Thus ip n (x) is recessive both at +00 and —00 and is unique up to 
a multiplicative constant. We can suppose that ho is sufficiently small so that il} n (x) 
has at least one zero in the interval By Theorem [TJ ip n (x) has no zeroes in 

the set R\ ]x-, x + [. 

Let us choose a point X\ <E]x_,x + [ independently of n. Let x[ be the zero of ip n 
which is nearest to x±. This means that x[ depends on n but the distance between x\ 
and x[ tends to zero as n tends to infinity. Denote by m + and m_ the number of zeroes 
of ip n in the interval [ [ and 1 ], respectively (hence the zero x[ is counted 
both in m + and m_). Denote by (+{x) the function defined by relations ( 11 lap and 
(lllbp . with x being replaced by x + . In virtue of Theorem [3, there exists a constant 
c + > (independent of n) such that 

IC+(X) -^n /3 am+| < 

m + 

for all n > uq. An application of the mean value theorem, 

|w 3//2 — f 3 / 2 | < - (max{«, w}) 1 ^ 2 \u — v\ for w > 0, v > 0, 
yields the inequality 

||C + K)| 3 / 2 -^|a m J 3 / 2 | <^Q£ + v^T^dxy /3 ^ (22) 

which is valid for all sufficiently large n. Analogously, for the other turning point we 
get the estimate 

|K-K)| 3/2 -^|a m _| 3 / 2 |<^(^ f X+ y^VixjdxT -^ (23) 

where again c_ > is a constant independent of n. Set 

/3 rx+ \i/3 

c = I — J \j~E — V{x) dx) max{c_,c + }. 
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Combining (1221) and (I23p we arrive at the inequality 



h„ 



E-V(x) dx-- (|< 



1 3/2 



+ \a r 



1 3/2 



< C 



m 



1/3 



+ 



m 



1/3 



Let m = m(n) be the number of zeroes of if> n (x). Obviously, m = m_ + m + — 1. 
Recalling the asymptotic behavior of h n , as stated in ( J20l) (see also (JTHl) ). as well as 
the asymptotic formulas (12~T1) for the roots of the Airy functions we finally find that 



n 



m 



0(n 



< 



71 



m 



1/3 



m 



1/3 



By Theorem [7J both m_ and m + tend to infinity as n tends to infinity. This implies 
that m(n) = n for all sufficiently large n and therefore, in virtue of Theorem [21 E 
is the nth eigenvalue of the Hamiltonian H(h n ) (with the numbering starting from 
n = 0). 

All estimates can be carried out in a uniform manner for E being replaced by A 
running over some neighborhood of E. We conclude that 

with the assumptions on V(x) formulated in the Introduction, there exist n £ N 
and a neighborhood Ue of E such that for all n > n and A £ Ue, A equals exactly the 
nth eigenvalue of H(h n (X)) (with h n (X) introduced in Theorem^). 



6 Proof of the formula 

Here we prove the limit (EI). We know that there exists a sequence of positive numbers, 
{^•n}5£L no > sucn ^ na ^ E ^ S ^ ne nth. eigenvalue of H{h n ) (Theorem E]). This sequence is 
strictly decreasing and tends to 0. We even known that H n ~ n^ 1 as n — > cxd (see 
( !20l) ). Therefore everywhere in what follows the symbol 0(h) should be understood as 
a substitute for 0(n _1 ). 

Let us fix xi, x[,x'( £ x+[, x[ < x\ < x". For a given h = h n we shall denote by 
^ a conveniently normalized eigenfunction corresponding to the eigenvalue E = E n . 
Hence ip is recessive both at +oo and — oo. The normalization is fixed by requiring 
the eigenfunction ip to coincide on the interval }x[, +oo[ with the solution described in 
Theorem [5] (with f(x) = V(x) — E and Xq = x + being the single turning point in this 
interval). Theorem [5] is also applicable to the interval ] — oo, x'{[ containing the turning 
point a;_. On this interval, ip equals k times the solution described in Theorem [5] for 
some K £ C \ {0}. 

There exists a neighborhood of E, Ue = ]E_, E + [, such that any A £ Ue satisfies 
the same assumptions as those imposed on E. Recall that we have fixed k £ Z. For all 
sufficiently large n, the (n+k)th eigenvalue of H(h n ), called E n+ k, exists and lies in Ue- 
For brevity we shall denote E n+ k sometimes by E. We show below that E — E = 0(h), 
see fflM . The eigenfunction of H(h n ) corresponding to the eigenvalue E = E n+k and 
coinciding on Jx^+oof with the solution from Theorem E] will be denoted by if). In 
this case, too, there exists k £ C \ {0} such that on the interval ] — oo, x'[ [ , if) equals 
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R times the solution from Theorem El Furthermore, denote by x± the turning points 
corresponding to E, i.e., V(x±) = E. Since V(x±) — V(x±) = E — E and V\x±) ^ 
it is clear that x± — x± = 0(h) as well. 

The verification of (jfJJ) is based on a series of estimates relying on Theorem [51 This 
will be done in several steps. 

(1) Relation between E and E. Let E m (h) be the mth eigenvalue of H(h). From the 
perturbation theory [TT] one deduces that if it exists and lies below Voo then E m (h) 
is strictly increasing and real analytic as a function of H. According to the conclusion 
of Section EJ E m {h) and h m (X) are mutually inverse functions. Therefore if h = h n (E) 
then h = h n+ k(E). Thus we have 

h n {E) = h n+k {E) 

and from the asymptotic formula (|2*U|) we get 

(2n + 2k + l)J(E) - (2n + 1) J(S) = 0(n -1 ). 

Since 

J(£) = J(E) + ^^(E - E) + 0((E' - E) 2 ) 



we finally arrive at the equation 

2k J(E) 
2n + 1 T(E) ' 

whose solution satisfies 

E = 



E + E = 0{n- 2 ) + 0((E - Ef) 



^ J(E) k 2 , 
E + -^-Jr - + 0(n~ 2 " 



- nr] (24) 

1 [E) n 

(2) Asymptotic behavior of k and k. On the interval Jx^x^ one can compare the 
asymptotics of the solutions which are respectively recessive at +oo and — oo and infer 
this way the asymptotic behavior of k as h — > 0. For a moment we shall distinguish by 
a subscript the functions (± related to the turning points x± and defined respectively 
on the intervals [x[, +oo[ and ] — oo,x"}. Thus 



1 2/3 



\f(t)\dt 



IC- 



1 2/3 



l/(*)|dt 



and both C+/f an d C-/f are positive functions on their domains. We have 

1/4 



for x > x[, and 



i/j(x) 



ip(x) 



(Ai(h^ 3 C + )+e + (h,x)) 
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for x < x". Suppose that ]. Recalling that 

Ai(-,) = -^(cosQ^-^ + 0(,- 3 / 2 )) asz^+oo (25) 
and the error term estimates from Theorem [5] we arrive at the equality 

cosQ trW - + 0(h) = k (cosQ ^iC-i 3 / 2 - ^ + O(ft)) . 
Furthermore, in virtue of ( !20l) it holds 

| (K + | 2/3 + IC-| 2/3 ) = £ dt = (n + vr + O(ft). 
Combining the last two equalities we find that 

cosQ ^K+r - + 0(h) = k ((-IfcosQ r^C+l 3 / 2 - ^ + O(ft)) . 
For 7i sufficiently small it clearly exists x G [x^x'/] such that 

cosQ^|C + | 3 / 2 -^=l. 

It follows immediately that 

«=(-l) n + 0(fi). (26) 

Similarly, 

k = (_i)™+ fc + 0(ft). (27) 

^ T/ie leading asymptotic term on the interval ]x + — 5, oof. Fix 5 > sufficiently 
small (at least x\ < x + — 5). Let us show that 

poo px—-\-S 

/ 4> 2 dx = 5 1 / 2 0(h 1 / 3 ), / ^dx = 6^ 2 0(h^ 3 ). (28) 

J X+—S J— oo 

We shall verify only the first equality in ( |28i) . In view of ( 1261) and ( 1271) . the verification 
of the second one is analogous. 

Here and everywhere in what follows the symbol 0(h e ) should be interpreted prop- 
erly. It means that there exists a constant c > (independent of 5) and ho(5) > 
such that for all h, < h < h (5), it holds \0(h e )\ < ch £ . 

First let us estimate the contribution from the leading asymptotic term of ip. 
Applying the substitution x = £(h 2 / 3 z) we get the expression 

poo / /-\ 1 / 2 POO 

/ {)) Ai(rn) 2 dx = ^ 3 / I Ai(,) 2 d,. (29) 

Jx+s \ fj Jh-VH^+s) f(i\ h ' z )) 
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By the assumptions, there exist x 2 > x + and C\ > such that f(x) > C\ for x > x 2 - The 
function ((x)/f(x) is continuous on the interval [xi,^] and therefore it is majorized 
on this interval by a constant c 2 > 0. This also means that 

< TTcT~V\ ^ C2 for ^ CO 1 *)- 



This way we get the following upper bound on (|29l) . namely 



ft 2 / 3 / c 2 Ai(z) 2 dz + / - Ai(z) 2 dz 

< c 2 h 2/3 (Ai'(x) 2 - x Ai(x) 2 ) + o(h 4/3 ). 

Here we have used the knowledge of the primitive function 

J Ai(x) 2 dx = x Ai(x) 2 - Ai'(x) 2 . 



In addition to formula (1251) let us recall also the asymptotic behavior of the derivative 
of the Airy function, 



Af(-z) = ^ (sin Q z*' 2 - + 0{z~^ as z - +oo. 



(30) 



Since C(a; + — 5) = —('(y)$ for some ?/ G — 5, s + ] we find that for £ = h 2 ^ 3 ((x + — 5) 
it holds 

|fr 2 / 3 Ai'(a;) 2 | < const h 2 ' 3 (7r 2 / 3 5) 1/2 = const h^H 1 ' 2 

and 

|/i 2 / 3 x Ai(x) 2 | < const h 2 ' 3 h- 2 / 3 5 (/r 2 / 3 5)~ 1/2 = const tVH 1 ' 2 . 
We have shown that 

1/2 



J QJ Ai(/l- 2 / 3 C) 2 dx = S^Otf' 3 ) 



T/ie error term on t/ie interval ]x + — 5, oo[. Further let us write 

1/2 



It is known that 



^ 2 =(£) Ai(/r 2 / 3 2 + e 2 (M). 

Ai(x) < —^=x~ 1/4 exp( -- h~ l x 312 J for x > 0, 
2^71 \ 3 / 

see [TBI Chapter 1 1] . Using also the estimates of error terms from Theorem one can 
check that 



\e 2 (h,x)\ < const /" 1/2 exp (-^ /T 1 ^ 72 ^ fr 4/3 /or x > 



x + . 
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It follows that 



e2(h, x) dx 



< const I™ f- 1 ' 2 exp ( h V ; ) d.r 
= const J™ exp (-1 h~ l y^ dy. (31) 

There exists c > such that for y > 0, f(£(y)) 1 < c(l + y _1 ). Hence (13~T1) is majorized 

by 

const ^(Z/ 172 + Z/ _1/2 ) exp f-~ rV 3/2 ) = 0(/l 5 / 3 ). 

The asymptotic formula ( 1251) implies that | Ai(x)| < const |x| _1//4 for x < 0. Re- 
calling once more Theorem [5] we have 



E2(h, X) dx 



J-i 



< const ft 4 / 3 / I/I" 1 / 2 dx = 0(h 4/3 ). 



(32) 



This concludes the verification of 
(5) Oscillating integral on the interval ]x±, x + — 5[. By the usual integration by parts 
one can verify the following claim. 

Let [a,b] be a compact interval, F G C l (\a,b]), fi G C 2 ([a,b]) and v{h,z) be twice 
continuously differentiable in z on [a,b]. Assume that //(z) nowhere vanishes on [a,b] 
and 

sup \d x v(h,z)\ = 0(1), sup \d 2 z v(h,z)\ = 0(1). 
Then for all sufficiently small H it holds true that 



f F(z) sinfi" 1 /J,(z) + u(H,z))dz 

J a 



< const h 



where the constant depends only on the length of the interval [a, b] and on the quantities 

ll^llc /^o 2 |l-^llc*||/ i "||c ) fJ'o ll^'llc 

with 

/io = min 

ze[a,b] 

and || ■ He* standing for the norm in the Banach space C([a, b}). 
As a consequence we find that if W G C l (1R) then 

_-<5 



w 



VE-V 



sin 



h-\\c\ 3/2 + \C\ 3/2 ))dx = 6- 1 0(h). 



(33) 



To show this asymptotics it suffices to set in the above claim F = W/ y/E — V, \i 
(4/3) (CI 3/2 and 



h 



-1 



x- + 



E — V(t) dt 
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'E - V(t) dt . 



Hence jx'(z) = -lyjE - V{z) and 



E-E 



d 2 u(H,z) = — — \^E-V(z) + \JE-V{z) 

d 2 Xh,z) = ^ V '(z)(E-V(z))- 1/2 (E-V(z))- 1/2 

x ^y/E-V(z) + yjE-V(z)^J . 

(6) The leading asymptotic term on the interval }x\,x + — 5[. Let us check the con- 
tribution to the matrix element coming from the interval [xi,x + — 5]. The leading 
asymptotic term in the expansion of if) is given in ( JT6l) . We also need the asymptotic 
behavior of the Airy function (1251) and the fact that the function //£ is continuous 
and hence bounded on the interval We conclude that 



if) r- 

Observe that 



£4/3 / = r l 0( £4/3 ); 



and on the interval [x±,x+ — 5], 

(e - vy l/i — (e — vy l/ \i + r 1 o(^)). 

From the boundedness of W and from an estimate similar to ( 1321) it follows that 

m; + -(5 ^x + -(5 / /■ \ 1/4 / > \ 1/<4 

y W^dx = y ^(jj (7) Ai(^ 2/3 C)Ai(^ 2/3 C)dx + 0(n 4/3 ) 

£1/3 W 

= —L w~ Al+6 0{h)) 

x cosQ - ^ oobQ h-^ 2 - |) da; 

+ r 1 o(^ 4/3 ). 

Using the asymptotic behavior ( 1331) we have 
-x+-5 _ ^1/3 / 2 



*W dx = ^ jT" cosQ ^(Kl 3 / 2 - |CI 3/2 )) dx + r 1 o(^ 3 ). 



(34) 

(7) The argument of the cosine on the interval ]x±,x + — 5[. Let us show that for 
x e [x\, x + — 5], 

I fr-WCl 3 ^ ~ \C\ 3/2 ) = ~ r(x) + 6^0(1) (35) 
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where t(x) was defined in (jSJ). We have 

^ _1 (lcr /a -ICI 3/2 ) = fr 1 (J* + VW^Vdt - £ VE-Vdt 

= ^{^JJ (y/E-V - VE-V^) dt 



+ / VE-Vdt- / VE-Vdt J. 



Set temporarily 



g(y) = I y/V(y)-V(t)dt. 



Then for y lying between x + and x + it holds 

W(y)\ = 4 / V ' (y) dt 



1 r y 

2 J x+ -s ^V(^V(t) 



1 /■» dt 

< - const / —== 
2 



< const \/|aJ+ — x+| + <5. 



Hence 



/ y/E^Vdt- VE-Vdt = \g(x+)-g(x 

Jx+—8 J xj^—S 



< const a/| 



,X_i_ — ,T_ 



— x_ 



Furthermore, 



V-E- V -VE-V 



(E — E) : 



E-E 

2 ^ E ~ V 2^Ve~^v + Ve-vJve~^v 

(E-E) 2 



< 



2(E - Vfl 2 



and 



From (j2U) it follows that 



tr\E-Et) = -£- + 0{h) 
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where T is the period of the classical motion, see (J7J). Altogether this means that 

/ 2ttA; A r+- s dt 



+ 6- 1/2 o(h) 

2nk r x + dt 



T J x 2^E-V{t) 
+ 5 1/2 0(l) + 5- 1/2 0(h) 

Relations flHol) and f[3"Tj) jointly imply 

(8) The final step. From ( !34"|) and (l35j) we derive that 

/ Wil)i/jdx — ' ' 



2^ - V(t) 

(37) 



/ .ri 



2tt U« V^-V(a;) 



X cos 



h l/3 / yx+ / 2?rfc 



COS 



r(z) j dx + 5 l/2 0{\) I . (38) 



2vr ^ ^E-y(x) V r 

The interval [x- + 5, x±] can be treated similarly. We have 

f X1 W^dx = KK—{ r W{X) cos('^T4x))dx + 5 1 / 2 0(l) 
J X .+5 2tt \J X _ ^ E - V(x) V T y >) 1 \ 

where 

T_(X) = - — = — — 1 — T[X). 

Taking into account also ( 1261) and ( 1271) we finally find that 

r +J ^ - ^ (r *>) + ■ w 

From the boundedness of W and relations ff28l . fl38|) and fl39|) it follows that 

dx = -— / W(g(t))e ifca; ' dt + 5 1/2 0(1) . (40) 



2vr 

As a particular case, with W(x) = 1 and k = 0, we have 

/" ^ dx = ^!( T + 5 i/2 0(1) ). (41) 
2tt 

The same relation holds also for the squared norm of ip. 

Relations ( 1401) and ( HIT) imply that there exists c > such that for all sufficiently 
small positive 5 and all n, n > no(5), it holds 

(n\W(x)\n + k) - i / ^(g(t))e ifc ^dt < c5 1/2 . 

Since 5 is arbitrary this concludes the verification of the limit ([6]). 
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